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Abstract 

This talk provides a basic introduction for students interested in the re- 
^ ' sponses of detectors to solar, supernova, and other low-energy neutrino sources. 

o 



Some of this nuclear physics is then applied in a discussion of nucleosynthesis 
within a Type II supernova, including the r-process and the //-process. 

O ' 

Q\ . 1. Introduction 

It is a pleasure to have this opportunity to visit Tokyo Metropolitan Uni- 
versity and address this group of students and researchers interested in neutrino 
physics. Professor Minakata has asked me to provide a pedagogical overview of 
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the nuclear physics governing the detection of solar, supernova, and other low- 
energy neutrinos. As the following presentation is very elementary, I apologize to 



$h ' those of you who are already familiar with the subject. 

The talk begins with a discussion of the allowed and first-forbidden re- 
sponses of nuclei to low-energy neutrinos. To illustrate how the allowed response 
can be crucial to efforts to detect solar neutrinos, I discuss the classic example 
of the 37 CI experiment. Similarly, first-forbidden responses are generally quite 
important to the interaction of heavy-flavor neutrinos from core-collapse super- 
novae. I discuss some examples from explosive nucleosynthesis - the r-process 
and the ^-process - to illustrate some of the issues. 



2. The allowed response 

Figure 1 shows several semileptonic weak interactions that take place be- 
tween nucleons or in nuclei [19]. Among such reactions important to astrophysics, 
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two of the most familiar are the decay of the free neutron 

n — > p + e~ + v e , 

a reaction that influences the n/p ratio in big-bang nucleosynthesis, and the driv- 
ing reaction of the solar pp chain 

p + p — > d + e + + v e . 

The latter can be thought of as the decay of a free proton in the plasma, made 
possible energetically by the proximity of a second proton, within the range of 
the nuclear force (several fermis), so that the final n+p state can form a bound 
deuteron. It is the binding energy of the deuteron that allows the reaction to take 
place. 




Fig. 1. Semileptonic weak interactions of interest: a) charged current neutrino 
reaction (e.g., 37 Cl(z/, e~) 37 Ar); b) f3~ decay; and c) electron capture (e.g., 
37 Ar(EC) 37 Cl). 



As preparation for our discussion of nuclei, first consider the rate for neu- 
tron (3 decay 

j \n/r\2 d^Pp M p (ftp e 7Tl e (ftpv VTL V 4 4 

dU=m (2^^(2^^(2^^ (27r) 6 ^-^-Pe-P,). (1) 

(Note that I use a spinor normalization convention where all fermions are treated 
as massive, including neutrinos.) The invariant amplitude M is taken to be a 
contact current-current interaction, because the momentum transfered between 
the leptons and nucleon is so much smaller than the mass of the W boson. Thus 

G 

M = cos 9 c -^u(p)Y (I - g A K)u(n)u(e)^(l - -y 5 )v(y) (2) 



3 



where Gp is the weak coupling constant measured in muon decay and cos^c 
gives the amplitude for the weak interaction to connect the u quark to its first- 
generation partner, the d quark. The origin of this effective amplitude is the 
underlying standard model predictions for the elementary quark and lepton cur- 
rents. The weak interactions at this level are predicted by the standard model to 
be exactly left handed. Experiment shows that the effective coupling of the W 
boson to the nucleon is governed by 

(l-^Ts) 

where qa ~ 1.26. The axial coupling is thus shifted from its underlying value 
by the strong interactions responsible for the binding of the quarks within the 
nucleon. 

The extension to nuclear systems traditionally begins with the observation 
that nucleons in the nucleus are rather nonrelativistic, v/c ~ 0.1. The (5 decay 
amplitude u{p)^{l — gAlb) u { n ) can be expanded in powers of p/M. The leading 
vector and axial operators are readily found to be 

Thus it is the time-like part of the vector current and the space-like part of the 
axial- vector current that survive in the nonrelativistic limit. 

(In a nucleus these currents must be corrected for the presence of meson 
exchange contributions. The corrections to the vector charge and axial three- 
current, which we just pointed out survive in the nonrelativistic limit, are of order 
(v/c) 2 ~ 1%. Thus the naive one-body currents are a very good approximation to 
the nuclear currents. In contrast, exchange current corrections to the axial charge 
and vector three-current operators are of order v/c, and thus of relative order 1. 
This difficulty for the vector three-current can be largely circumvented, because 
current conservation as embodied in the generalized Siegert's theorem allows one 
to rewrite important parts of this operator in terms of the vector charge operator. 
In the long-wavelength limit appropriate to f3 decay, all terms unconstrained by 
current conservation do not survive. In effect, one has replaced a current operator 
with large two-body corrections by a charge operator with only small corrections. 
In contrast, the axial charge operator is significantly altered by exchange currents 
even for long-wavelength processes like (3 decay. Typical axial-charge f3 decay 
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rates are enhanced by ~ 2 because of exchange currents.) 

If such a nonrelativistic reduction is done for our nucleon (3 decay ampli- 
tude, one obtains 

G 

M ~ cos6» c -^ (0 t (p)0(n)M(e)7°(l - 7 5 )u(i/) - $(p)g A v(f>(n) ■ u(e)^(l - ^)v{yfj 

(3) 

where the are now two-component Pauli spinors for the nucleons. The above 
result is written for the /3 decay n — > p. It is convenient to generalize it for p <-> n 
by introducing the isospin operators r± where r + | n) = | p) and r_| p) = | n), 
with all other matrix elements being zero: the free proton does not (3 decay, of 
course, but this is good preparation for the generalization to nuclei. Finally, we 
square the invariant amplitude, integrate over the outgoing electron, neutrino, 
and final nucleon three-momenta, average over initial nucleon spin, and sum over 
final nucleon spin, electron spin, and neutrino spin. The result is 

W = G#cos 2 6> c ^ jj w - e ^ e2 - m " de \ (K/H^IIOI' + ^K/II^IIOI 2 ) (4) 

where / and i are the final and initial nucleon states, w is the energy release in 
the decay, and e is the electron energy. The r+ operator corresponds to f3~ decay 
and the r_ to (3 + decay. The notation || denotes a matrix element reduced in 
angular momentum. One immediately sees, for large energy release w, that rates 
scale as w 5 . 

This result easily generalizes to nuclear decay. Given our comments about 
exchange currents, the first step is the replacement 

A 

1=1 

A 

i=i 

We also have to worry about an approximation in our nucleon (3 decay discussion, 
the treatment of the nucleon as an elementary, structureless particle. This is 
certainly appropriate for momentum scales below the inverse size of the nucleon, 
as the nucleon's structure then cannot be resolved, and for energy transfers small 
compared to nucleon excitation energies. Both conditions are easily satisfied in 
neutron (3 decay. In nuclear (3 decay the issue is not so clear, especially as decays 
can often populate a collection of states in the daughter nucleus. If the lepton 
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states are treated as plane waves, the operators further generalize to 

A 

i=i 
i=i 

where f(i) is the coordinate of the ith nucleon relative to the nuclear center of 
mass. (The center-of-mass coordinate would be integrated out to give the overall 
three-momentum conservation for the f3 decay.) In f3 decay and in solar neutrino 

— * 

reactions, the three-momentum transfer to the nucleus \k\ is much smaller than 
the typical inverse nuclear size, ~ 160/A 1 / 3 MeV. Thus as long as one is interested 
in transitions where the operators r or ar connect the initial and final states of 
interest, the effects of the momentum transfer can be ignored. Of course, if this is 
not the case, then the transition amplitude is nonzero only because of the finite 
momentum transfer. If one expands the plane wave in powers of k ■ r*(i), then 
a transition has a degree of "forbiddenness" according to the number of powers 
required to produce a nonzero amplitude. 

For the moment we will restrict ourselves to allowed transitions where the 
effects of the momentum transfer can be ignored. The nuclear decay rate is then 
obtained by substituting into the neutron result 

A A 

^t-t(I^IIE^WIK)I 2 + ^I(/IIE^K«IK)I 2 )- (5) 

ZJ i 1 i=l i=l 

The factor 1/(2 Jj + 1) replacing the 1/2 in the neutron result comes from the 
average over initial nuclear spin directions. As the nuclear Coulomb field can 
significantly distort the wave function of the outgoing electron or positron, a final 
step is to correct the lepton phase space by 

F(Z,e) = \F (Z,e)f = J^- i where V = 

where (3 is the electron/positron velocity and F (Z, e) is the s-wave Coulomb 
wave function in the field of the daughter nucleus of charge Zf, evaluated at 
the nuclear origin. (This is a reasonable approximation for small Zf] for heavier 
nuclei, however, the usual procedure is to solve the Dirac equation for an exten- 
sive nuclear charge, evaluating the resulting wave function at the nuclear surface.) 



The spin-independent and spin-dependent operators appearing above are 
known as the Fermi and Gamow- Teller operators. The Fermi operator is propor- 
tional to the isospin raising/lowering operator: in the limit of good isopsin, which 
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typically is good to 5% or better in the description of low-lying nuclear states, 
it can only connect states in the same isospin multiplet, that is, states with a 
common spin-spatial structure. If the initial state has isospin (T«, M^), this final 
state has (T,, Mti ± 1) for f3~ and j3 + decay, respectively, and is called the isospin 
analog state (IAS). In the limit of good isospin the sum rule for this operator in 
then particularly simple 

The excitation energy of the IAS relative to the parent ground state can be 
estimated accurately from the Coulomb energy difference [9] 

1 728^ 

E ™ ~ W/» + 0.78 - L293)MeV ' (7) 

The angular distribution of the outgoing electron for a pure Fermi (N, Z) + v — > 
(N — 1, Z + 1) + e~ transition is 1 + f3 cos9 ue , and thus forward peaked. Here (3 
is the electron velocity. 

The Gamow- Teller (GT) response is more complicated, as the operator can 
connect the ground state to many states in the final nucleus. In general we do 
not have a precise probe of the nuclear GT response apart from weak interactions 
themselves. However a good approximate probe is provided by forward-angle 
(p,n) scattering off nuclei, a technique that has been developed in particular by 
experimentalists at the Indiana University Cyclotron Facility. The (p,n) reaction 
transfers isospin and thus is superficially like (i/,e~). At forward angles (p,n) 
reactions involve negligible three-momentum transfers to the nucleus. Thus the 
nucleus should not be radially excited. It thus seems quite plausible that forward- 
angle (p,n) reactions probe the isospin and spin of the nucleus, the macroscopic 
quantum numbers, and thus the Fermi and GT responses. For typical transitions, 
the correspondence between (p,n) and the weak GT operators is believed to be 
accurate to about 10%. Of course, in a specific transition, much larger discrep- 
ancies can arise. 

The (p,n) studies demonstrate that the GT strength tends to concentrate 
in a broad resonance centered at a position 5 = E GT — Ejas relative to the IAS 
given by [14] 

N - Z 

5 ~ (7.0 - 28.9 — — ) MeV. (8) 

Thus while the peak of the GT resonance is substantially above the IAS for N ~ Z 
nuclei, it drops with increasing neutron excess. Thus 5 ~ for Pb. A typical 
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value for the full width at half maximum T is ~ 5 MeV. 

The approximate Ikeda sum rule constrains the difference in the f3~ and 
j3 + strengths 

E(I^?t(/T)| 2 - \Mg T ({3 + )\ 2 ) = 3(iV - Z) (9) 
/ 

where 

iM^ni^yyVrK/llE^KWIN)! 2 - (io) 
zj « ' 1 i=i 

In many cases of interest in heavy nuclei, the strength in the f3 + direction is 
largely blocked. For example, in a naive 2sld shell model description of 37 C1, 
discussed below, the p — > n direction is blocked by the closed neutron shell at 
N=20. Thus this relation can provide an estimate of the total f3~ strength. Ex- 
periment shows that the f3~ strength found in and below the GT resonance does 
not saturate the Ikeda sum rule, typically accounting for ~ (60 — 70) % of the 
total. Measured and shell model predictions of individual GT transition strengths 
tend to differ systematically by about the same factor. Presumably the missing 
strength is spread over a broad interval of energies above the GT resonance. This 
is not unexpected if one keeps in mind that the shell model is an approximate 
effective theory designed to describe the long wavelength modes of nuclei: such 
a model should require effective operators, renormalized from their bare values. 
Phenomenologically, the shell model seems to require g e / ~ 1.0 as well as a small 
spin-tensor term (a <g> Y 2 (f))j = i of relative strength ~ 0.1 [2]. 

The angular distribution of GT (TV, Z) + is e -> (N - 1, Z + l) + e~ reactions 
is 3 — j3cos9 ue , corresponding to a gentle peaking in the backward direction. 

The above discussion of allowed responses can be repeated for neutral cur- 
rent processes such as (u, v'). The analog of the Fermi operator contributes only 
to elastic processes, where the standard model nuclear weak charge is approxi- 
mately the neutron number. As this operator does not generate transitions, it is 
not yet of much interest for solar or supernova neutrino detection, though there 
are efforts to develop low-threshold detectors (e.g., cryogenic technologies) where 
the modest recoil nuclear energies might be detectable. The analog of the GT 
response involves 

\MgrW)\ 2 = ^VtK/IIE^)^!^)! 2 - (") 

The operator appearing in this expression is familiar from magnetic moments and 
magnetic transitions, where the large isovector magnetic moment (//„ ~ 4.706) 
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often leads to it dominating the orbital and isoscalar spin operators. 

3. The Response of the 37 C1 Detector 

An interesting example of these issues in the context of a practical detector 
is provided by the 37 C1 solar neutrino experiment of Davis and his collaborators. 
Davis succeeded in recovering and counting the few atoms of 37 Ar produced by 
solar neutrinos in a 0.615 kiloton C2CI4 detector via the reaction 37 Cl(z/, e~) 37 Ar. 
The capture rate determined from three decades of measurements in 2.56 ± 0.16 
± 0.16 SNU [7] (1 SNU = 10~ 36 captures/target atom/s), or about 1/3 that 
predicted by the standard solar model, conventional particle physics, and the ac- 
cepted value for the 37 C1 neutrino capture cross section. This experiment was the 
first manifestation of the solar neutrino problem and remains crucial to current 
conclusions that neutrino oscillations may be responsible for the neutrino deficit. 

The strong conclusions drawn from the 37 C1 experiment depend on an ac- 
curately determined neutrino capture cross section. Because the threshold for 
37 Cl(z/, e _ ) 37 Ar (0.814 MeV) is well above the pp neutrino endpoint, the impor- 
tant neutrino sources are from the 7 Be and 8 B solar reactions. The 7 Be neutrinos 
can only excite the transition to the ground state of 37 Ar, which is relatively weak 
(log/t = 5.10). Thus the capture rate should be dominated by the high energy 
8 B neutrinos (endpoint ~ 15 MeV). 

The nuclear (not atomic) mass difference between 37 CI and 37 Ar is 0.303 
MeV. The Coulomb energy difference formula (Eq. 7) for the position of the IAS 
gives 

E IAS ~ 5.22MeV. 

So we conclude that the analog state should reside at ~ 4.92 MeV in 37 Ar. Ex- 
periment has identified the IAS at 4.99 MeV. In the limit of good isospin the 
superallowed (Fermi) transition to the IAS has |M^| 2 = N - Z = 3.0; this transi- 
tion accounts for about 70% of the 8 B capture rate. 

Now the interesting issue is the model-dependent GT response. While we 
have noted that the total GT response is about three times the Fermi response 
(taking g e /^ ~ 1), its contribution to the capture rate depends on its distribution, 
particularly at low excitation energies where the 8 B neutrino cross section phase 
space is large. As the effective particle breakup threshold for 37 Ar is 8.79 MeV, 
GT transitions to states above this energy clearly do not contribute. According 
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to our estimate (see Eq. (8)) for 5 = E GT — E IAS ~ 4.66 MeV, the peak of the 
GT distribution should be at an excitation energy ~ 9.6 MeV, relative to the 
ground state of 37 Ar. The two strongest peaks in the forward-angle (p,n) studies 
are in the region between 7 and 10 MeV, roughly in accord with expectations. 
Thus much of the GT strength is in the continuum, and still more resides above 
the IAS, where the neutrino phase space drops rapidly with increasing excitation 
energy. 

To put the capture cross section on firm ground, a reliable map of the 
strength and distribution of the GT bound state response is needed. In 1964 
Bahcall and Barnes [4] pointed out that the needed information could be obtained 
from the delayed proton spectrum following the j3 decay of 37 Ca, as illustrated in 
Fig. 2. Assuming isospin invariance, the decay 37 Ca(/3 + ) 37 K is the mirror reaction 
to 37 Cl(z/, e _ ) 37 Ar. As the 37 K levels above the first excited state are unstable to 
proton emission, the allowed matrix elements for these levels can be deduced from 
the spectrum and intensities of the delayed protons. The transition to the ground 
state of Ar is known directly, as this transition determines the electron capture 
lifetime of 37 Ar. The final needed constraint on the transition to the first excited 
state is imposed by the total rate for 37 Ca (3 decay. Thus, to the extent that 
isospin invariance relates the mirror systems accurately, the needed GT strengths 
can be taken entirely from experiment. 

The 37 Ca(/9 + ) 37 K delayed proton spectrum was measured [16]] by two 
groups; the deduced ft values were the basis for the 37 C1 cross section used for 20 
years. Interestingly these early experiments were flawed because of a simplifying 
assumption, that the delayed protons from 37 K were accompanied by production 
of the daughter nucleus 36 Ar in its ground state. In 1987 Adelberger and Haxton 
[1], noticing that the GT distribution deduced from the delayed proton experi- 
ments differed significantly from that recently measured [17] in 37 Cl(p,n), argued 
that the likely source of this discrepancy was the population of 36 Ar in its 2 + first 
excited state (1.97 MeV) in the delayed proton experiments. The states popu- 
lated in 37 K by allowed (3 decay have the spins and parity l/2 + , 3/2 + , and 5/2 + . 
Thus the reason that the 36 Ar first excited state should be important is clear: the 
3/2 + and 5/2 + states can populate the 2 + state by s-wave proton emission, while 
the ground state requires d-wave emission. 

The conclusion was that the 37 Ca experiment had to be redone in a kine- 
matically complete way, where 1.97 MeV 7s accompanying the decay of the 2 + 
state could be observed in coincidence with the delayed protons. A series of 
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Fig. 2. Decay schemes of A=37 nuclei. For simplicity only the lowest T = 1/2 
and T = 3/2 levels are shown. The important particle decay thresholds for 
states in 37 Ar and 37 K are indicated. Atomic masses in MeV are indicated by 
[0.00], etc. All other energies are excitation energies in the indicated nuclei. 

elegant experiments were conducted by Garcia et al.[10], resulting in the determi- 
nation (t( 8 B) = 1.09 ± 0.09. Interestingly, this value was little changed from that 
used previously: ignoring the population of the 2 + state produced two largely 
compensating errors. The affected transitions were placed too low in energy (by 
1.97 MeV), but their strengths were also underestimated as the wrong 37 Ca (3 
decay phase space was then employed. However, the sizeable discrepancies be- 
tween the 37 Ca /3 decay and (p,n) mappings were largely resolved, thus restoring 
confidence that the capture rate uncertainties in the Davis experiment were under 
control. 

The reason for the discussions of this section is to illustrate that a reli- 
able cross section for the 37 C1 experiment was obtained only after complementary 
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calibration techniques were proposed, developed, and cross checked. This careful 
nuclear physics is a cornerstone of today's arguments that the solar neutrino puz- 
zle is likely due to new neutrino phenomena. 

4. Supernovae and Supernova Neutrinos 

Consider a massive star, in excess of 10 solar masses, burning the hydrogen 
in its core under the conditions of hydrostatic equilibrium. When the hydrogen is 
exhausted, the core contracts until the density and temperature are reached where 
3a -^ l2 C can take place. The He is then burned to exhaustion. This pattern (fuel 
exhaustion, contraction, and ignition of the ashes of the previous burning cycle) 
repeats several times, leading finally to the explosive burning of 28 Si to Fe. For a 
heavy star, the evolution is rapid: the star has to work harder to maintain itself 
against its own gravity, and therefore consumes its fuel faster. A 25 solar mass 
star would go through all of these cycles in about 7 My, with the final explosive 
Si burning stage taking a few days. The result is an "onion skin" structure of the 
precollapse star in which the star's history can be read by looking at the surface 
inward: there are concentric shells of H, 4 He, 12 C, 16 and 20 Ne, 28 Si, and 56 Fe at 
the center. 

The source of energy for this evolution is nuclear binding energy. A plot 
of the nuclear binding energy 5 as a function of nuclear mass shows that the 
minimum is achieved at Fe. In a scale where the 12 C mass is picked as zero: 



Thus once the Si burns to produce Fe, there is no further source of nuclear energy 
adequate to support the star. So as the last remnants of nuclear burning take 
place, the core is largely supported by degeneracy pressure, with the energy gen- 
eration rate in the core being less than the stellar luminosity. The core density is 
about 2 xlO 9 g/cc and the temperature is kT ~ 0.5 MeV. 



12 C 
16Q 

28 Si 
40 Ca 

56p e 

72 Ge 
98 Mo 



S I nucleon 
5 1 nucleon 
5/nucleon : 
S I nucleon 
5 1 nucleon 
S I nucleon 
S I nucleon 



0.000 MeV 
-0.296 MeV 
-0.768 MeV 
-0.871 MeV 
-1.082 MeV 
-1.008 MeV 
: -0.899 Mev 



Thus the collapse that begins with the end of Si burning is not halted by 
a new burning stage, but continues. As gravity does work on the matter, the 
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collapse leads to a rapid heating and compression of the matter. As the nucleons 
in Fe are bound by about 8 MeV, sufficient heating can release as and a few 
nucleons. At the same time, the electron chemical potential is increasing. This 
makes electron capture on nuclei and any free protons favorable, 

e~ + p — > v e + n. 

Note that the chemical equilibrium condition is 

He + Hp = + (E u ). 

Thus the fact that neutrinos are not trapped plus the rise in the electron Fermi 
surface as the density increases, lead to increased neutronization of the matter. 
The escaping neutrinos carry off energy and lepton number. Both the electron 
capture and the nuclear excitation and disassociation take energy out of the elec- 
tron gas, which is the star's only source of support. This means that the collapse 
is very rapid. Numerical simulations find that the iron core of the star (~ 1.2-1.5 
solar mases) collapses at about 0.6 of the free fall velocity [13]. 



In the early stages of the infall the u e s readily escape. But neutrinos are 
trapped when a density of ~ 10 12 g/cm 3 is reached. At this point the neutrinos be- 
gin to scatter off the matter through both charged current and coherent neutral 
current processes. The neutral current neutrino scattering off nuclei is partic- 
ularly important, as the scattering cross section is off the total nuclear weak 
charge, which is approximately the neutron number. This process transfers very 
little energy because the mass energy of the nucleus is so much greater than the 
typical energy of the neutrinos. But momentum is exchanged. Thus the neutrino 
"random walks" out of the star. When the neutrino mean free path becomes 
sufficiently short, the "trapping time" of the neutrino begins to exceed the time 
scale for the collapse to be completed. This occurs at a density of about 10 12 
g/cm 3 , or somewhat less than 1% of nuclear density. After this point, the energy 
released by further gravitational collapse and the star's remaining lepton number 
are trapped within the star. 

If we take a neutron star of 1.4 solar masses and a radius of 10 km, an 
estimate of its binding energy is 

CM 2 

— — ~ 2.5 x 10 53 ergs. 
2ir 

Thus this is roughly the trapped energy that will later be radiated in neutrinos. 
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The trapped lepton fraction Y L is a crucial parameter in the explosion 
physics: a higher trapped Y L leads to a larger homologous core, a stronger shock 
wave, and easier passage of the shock wave through the outer core, as will be dis- 
cussed below. Most of the lepton number loss of an infalling mass element occurs 
as it passes through a narrow range of densities just before trapping. The reasons 
for this are relatively simple: on dimensional grounds weak rates in a plasma go 
as T 5 , where T is the temperature. Thus the electron capture rapidly turns on 
as matter falls toward the trapping radius, and lepton number loss is maximal 
just prior to trapping. Inelastic neutrino reactions have an important effect on 
these losses, as the coherent trapping cross section goes as El and is thus least 
effective for the lowest energy neutrinos. As these neutrinos escape, inelastic reac- 
tions repopulate the low energy states, allowing the neutrino emission to continue. 

The velocity of sound in matter rises with increasing density. The inner 
homologous core, with a mass Mhc ~ 0.6 — 0.9 solar masses, is that part of the 
iron core where the sound velocity exceeds the infall velocity. This allows any 
pressure variations that may develop in the homologous core during infall to even 
out before the collapse is completed. As a result, the homologous core collapses as 
a unit, retaining its density profile. That is, if nothing were to happen to prevent 
it, the homologous core would collapse to a point. 

The collapse of the homologous core continues until nuclear densities are 
reached. As nuclear matter is rather incompressible (~ 200 MeV/f 3 ), the nuclear 
equation of state is effective in halting the collapse: maximum densities of 3-4 
times nuclear are reached, e.g., perhaps 6 • 10 14 g/cm 3 . The innermost shell of 
matter reaches this supernuclear density first, rebounds, sending a pressure wave 
out through the homologous core. This wave travels faster than the infalling 
matter, as the homologous core is characterized by a sound speed in excess of 
the infall speed. Subsequent shells follow. The resulting series of pressure waves 
collect near the sonic point (the edge of the homologous core). As this point 
reaches nuclear density and comes to rest, a shock wave breaks out and begins its 
traversal of the outer core. 

Initially the shock wave may carry an order of magnitude more energy than 
is needed to eject the mantle of the star (less than 10 51 ergs). But as the shock 
wave travels through the outer iron core, it heats and melts the iron that crosses 
the shock front, at a loss of ~ 8 MeV/nucleon. The enhanced electron capture 
that occurs off the free protons left in the wake of the shock, coupled with the 
sudden reduction of the neutrino opacity of the matter (recall a C oherent ~ N 2 ), 
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greatly accelerates neutrino emission. This is another energy loss. [Many nu- 
merical models predict a strong "breakout" burst of u e s in the few milliseconds 
required for the shock wave to travel from the edge of the homologous core to 
the neutrinosphere at p ~ 10 12 g/cm 3 and r ~ 50 km. The neutrinosphere is 
the term from the neutrino trapping radius, or surface of last scattering.] The 
summed losses from shock wave heating and neutrino emission are comparable to 
the initial energy carried by the shock wave. Thus most numerical models fail to 
produce a successful "prompt" hydrodynamic explosion. 

Two explosion mechanisms were seriously considered in the last two decades. 
In the prompt mechanism [8] described above, the shock wave is sufficiently strong 
to survive the passage of the outer iron core with enough energy to blow off the 
mantle of the star. The most favorable results were achieved with smaller stars 
(less than 15 solar masses) where there is less overlying iron, and with soft equa- 
tions of state, which produce a more compact neutron star and thus lead to more 
energy release. In part because of the lepton number loss problems discussed ear- 
lier, now it is widely believed that this mechanism fails for all but unrealistically 
soft nuclear equations of state. 

The delayed mechanism [5] begins with a failed hydrodynamic explosion; 
after about 0.01 seconds the shock wave stalls at a radius of 200-300 km. It ex- 
ists in a sort of equilibrium, gaining energy from matter falling across the shock 
front, but loosing energy to the heating of that material. However, after perhaps 
0.5 seconds, the shock wave is revived due to neutrino heating of the nucleon 
"soup" left in the wake of the shock. This heating comes primarily from charged 
current reactions off the nucleons in that nucleon gas; quasielastic scattering also 
contributes. This high entropy radiation-dominated gas may reach two MeV in 
temperature. The pressure exerted by this gas helps to push the shock outward. 
It is important to note that there are limits to how effective this neutrino energy 
transfer can be: if matter is too far from the core, the coupling to neutrinos is too 
weak to deposite significant energy. If too close, the matter may be at a temper- 
ature (or soon reach a temperature) where neutrino emission cools the matter as 
fast or faster than neutrino absorption heats it. The term "gain radius" is used 
to describe the region where useful heating is done. 

This subject is still controversial and unclear. The problem is numerically 
challenging, forcing modelers to handle the difficult hydrodynamics of a shock 
wave; the complications of the nuclear equation of state at densities not yet ac- 
cessible to experiment; modeling in two or three dimensions; handling the slow 



15 



diffusion of neutrinos; etc. Not all of these aspects can be handled reasonably 
at the same time, even with existing supercomputers. Thus there is considerable 
disagreement about whether we have any supernova model that succeeds in eject- 
ing the mantle. 

However the explosion proceeds, there is agreement that 99% of the 3 -10 53 
ergs released in the collapse is radiated in neutrinos of all flavors. The time scale 
over which the trapped neutrinos leak out of the protoneutron star is about 3 
seconds. (Fits to SN1987A give, assuming an exponential cooling e~^ T , r ~ 4.5 s 
[3]) Through most of their migration out of the protoneutron star, the neutrinos 
are in flavor equilibrium 

e.g., v e + i/ e <-+ i/p + Dp. 

As a result, there is an approximate equipartition of energy among the neutrino 
flavors. After weak decoupling, the z/ e s and z7 e s remain in equilibrium with the 
matter for a longer period than their heavy-flavor counterparts, due to the larger 
cross sections for scattering off electrons and because of the charge-current reac- 
tions 

u e + n <-> p + e~ 
u e + p <-> n + e + . 

Thus the heavy flavor neutrinos decouple from deeper within the star, where 
temperatures are higher. Typical calculations yield 

T Vli ~ T„ T ~ 8MeV 

T Ue ~ 3.5MeV ~ 4.5MeV. 

The difference between the u e and z/ e temperatures is a result of the neutron rich- 
ness of the matter, which enhances the rate for charge- current reactions of the 
u e s, thereby keeping them coupled to the matter somewhat longer. 

This temperature hierarchy is crucially important to nucleosynthesis and 
also to possible neutrino oscillation scenarios. The three-flavor MSW level-crossing 
diagram is shown in Fig. 3. One very popular scenario attributes the solar neu- 
trino problem to <-> v e transmutation; this means that a second crossing with 
a v T could occur at higher density. It turns out plausible seasaw mass patterns 
suggest a u T mass on the order of a few eV, which would be interesting cosmolog- 
ically. The second crossing would then occur outside the neutrino sphere, that is, 
after the neutrinos have decoupled and have fixed spectra with the temperatures 
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given above. Thus a v e <-> v T oscillation would produce a distinctive T ~ 8 MeV 
spectrum of u e s. This has dramatic consequences for terrestrial detection and for 
nucleosynthesis in the supernova. 




Fig. 3. Three-flavor neutrino level-crossing diagram. One popular scenario as- 
sociates the solar neutrino problem with u e <-> oscillations and predicts a 
cosmologically interested massive v T with v e <-> u T oscillations near the super- 
nova neutrinosphere. 



5. First Forbidden Responses and the Neutrino Process 

Core-collapse supernovae are one of the major engines driving galactic 
chemical evolution, producing and ejecting the metals that enrich our galaxy. 
The discussion of the previous section described the hydrostatic evolution of a 
presupernova star in which large quantities of the most abundant metals (C, O, 
Ne, ...) are synthesized and later ejected during the explosion. During the pas- 
sage of the shock wave through the star's mantle, temperature of ~ (1 — 3) ■ 10 9 K 
and are reached in the silicon, oxygen, and neon shells. This shock wave heating 
induces (7,01) <-> (01,7) and related reactions that generate a mass flow toward 
highly bound nuclei, resulting in the synthesis of iron peak elements as well as 
less abundant odd-A species. Rapid neutron-induced reactions are thought to 
take place in the high-entropy atmosphere just above the mass cut, producing 
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about half of the heavy elements above A ~ 80. This is the subject of the next 
section. Finally, the //-process described below is responsible for the synthesis of 
rare species such as n B and 19 F. This process involves the response of nuclei at 
momentum transfers where the allowed approximation is no longer valid. Thus 
we will use the //-process in this section to illustrate some of the relevant nuclear 
physics. 

One of the problems - still controversial - that may be connected with the 
neutrino process is the origin of the light elements Be, B and Li, elements which 
are not produced in sufficient amounts in the big bang or in any of the stellar 
mechanisms we have discussed. The traditional explanation has been cosmic ray 
spallation interactions with C, O, and N in the interstellar medium. In this pic- 
ture, cosmic ray protons collide with C at relatively high energy, knocking the 
nucleus apart. So in the debris one can find nuclei like 10 B, n B, and 7 Li. 

But there are some problems with this picture. First of all, this is an 
example of a secondary mechanism: the interstellar medium must be enriched in 
the C, O, and N to provide the targets for these reactions. Thus cosmic ray spal- 
lation must become more effective as the galaxy ages. The abundance of boron, 
for example, would tend to grow quadratically with metalicity, since the rate of 
production goes linearly with metalicity. But observations, especially recent mea- 
surements with the Hubble Space Telescope (HST), find a linear growth in the 
boron abundance [18]. 

A second problem is that the spectrum of cosmic ray protons peaks near 1 
GeV, leading to roughly comparable production of the two isotopes 10 B and n B. 
That is, while it takes more energy to knock two nucleons out of carbon than one, 
this difference is not significant compared to typical cosmic ray energies. More 
careful studies lead to the expectation that the abundance ratio of n B to 10 B 
might be ~ 2. In nature, it is greater than 4. 

Fans of cosmic ray spallation have offered solutions to these problems, e.g., 
similar reactions occurring in the atmospheres of nebulae involving lower energy 
cosmic rays. As this suggestion was originally stimulated by the observation of 
nuclear 7 rays from Orion, now retracted, some of the motivation for this scenario 
has evaporated. Here I focus on an alternative explanation, synthesis via neutrino 
spallation. 

Previously we spoke about weak interactions in nuclei involving the Gamow- 
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Teller (spin-flip) and Fermi operators. These are the appropriate operators when 
one probes the nucleus at a wavelength - that is, at a size scale - where the nu- 
cleus responds like an elementary particle. We can then characterize its response 
by its macroscopic quantum numbers, the spin and charge. On the other hand, 
the nucleus is a composite object and, therefore, if it is probed at shorter length 
scales, all kinds of interesting radial excitations will result, analogous to the vi- 
brations of a drumhead. For a reaction like neutrino scattering off a nucleus, the 
full operator involves the additional factor 

— * 

where the expression on the right is valid if the magnitude of k is not too large. 
Thus the full charge operator includes a "first forbidden" term 

A 
i=l 

and similarly for the spin operator 

A 

Y^fi <8>£(i)]j=o,i,2T3(i). 
i=i 

These operators generate collective radial excitations, leading to the so-called 
"giant resonance" excitations in nuclei. The giant resonances are typically at an 
excitation energy of 20-25 MeV in light nuclei. One important property is that 
these operators satisfy a sum rule (Thomas- Reiche-Kuhn) of the form 

A N7 A 

EK/IEKO^OIOI 2 -^-^ 

/ i=l 

where the sum extends over a complete set of final nuclear states. These first- 
forbidden operators tend to dominate the cross sections for scattering the high 
energy supernova neutrinos (i^s and v T s), with E u ~ 25 MeV, off light nuclei. 
From the sum rule above, it follows that nuclear cross sections per target nucleon 
are roughly constant. 

The El giant dipole mode described above is depicted qualitatively in Fig. 
4a. This description, which corresponds to an early model of the giant reso- 
nance response by Goldhaber and Teller [11], involves the harmonic oscillation 
of the proton and neutron fluids against one another. The restoring force for 
small displacements would be linear in the displacement and dependent on the 
nuclear symmetry energy. There is a natural extension of this model to weak in- 
teractions, where axial excitations occur. For example, one can envision a mode 
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similar to that of Fig. 4a where the spin-up neutrons and spin-down protons 
oscillate against spin-down neutrons and spin-up protons, the spin-isospin mode 
of Fig. 4b. This mode is one that arises in a simple SU(4) extension of the 
Goldhaber- Teller model, derived by assuming that the nuclear force is spin and 
isospin independent, at the same excitation energy as the El mode. In full, the 
Goldhaber- Teller model predicts a degenerate 15-dimensional supermultiplet of 
giant resonances, each obeying sum rules analogous to the TRK sum rule. While 
more sophisticated descriptions of the giant resonance region are available, of 
course, this crude picture is qualitatively accurate. 




Fig. 4. Schematic illustration of a) the El giant dipole mode familiar from 
electromagnetic interactions and b) a spin-isospin giant dipole mode associated 
with the first-forbidden weak axial response. 

This nuclear physics is important to the //-process [21]. The simplest ex- 
ample of //-process nucleosynthesis involves the Ne shell in a supernova. Because 
of the first-forbidden contributions, the cross section for inelastic neutrino scatter- 
ing to the giant resonances in Ne is ~ 3 ■ 10~ 41 cm 2 /flavor for the more energetic 
heavy-flavor neutrinos. This reaction 

v + A -> v' + A* 

transfers an energy typical of giant resonances, ~ 20 MeV. A supernova releases 
about 3 xlO 53 ergs in neutrinos, which converts to about 4 x 10 57 heavy flavor 
neutrinos. The Ne shell in a 20 M Q star has at a radius ~ 20,000 km. Thus the 
neutrino fluence through the Ne shell is 

4 • 1 57 

/ i , 9 ~ 10 38 /cm 2 . 

r 4vr(20, 000km) 2 1 

Thus folding the fluence and cross section, one concludes that approximately 
1 /300th of the Ne nuclei interact. 
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This is quite interesting since the astrophysical origin of 19 F had not been 
understood. The only stable isotope of fluorine, 19 F has an abundance 

19 F 1 
20 Ne ~ 3100' 

This leads to the conclusion that the fluorine found in toothpaste was created by 
neutral current neutrino reactions deep inside some ancient supernova. 

The calculation [21] of the final 19 F/ 20 Ne ratio is more complicated than 
the simple 1/300 ratio given above: 

• When Ne is excited by ~ 20 MeV through inelastic neutrino scattering, it breaks 
up in two ways 

20 Ne(z/, z/') 20 Ne* ^ 19 Ne + n ^ 19 F + e + + v c + n 
20 Ne(z/, z/) 2 °Ne* _>is p + p 

with the first reaction occurring half as frequently as the second. As both channels 
lead to 19 F, we have correctly estimated the instantaneous abundance ratio in the 
Ne shell of 

19 F _L_ 
20 Ne ~ 300' 

• We must also address the issue of whether the produced 19 F survives. In the 
first 10 -8 sec the coproduced neutrons in the first reaction react via 

15 0(n,p) 15 N 19 Ne(n,a) 16 20 Ne(n, 7 ) 21 Ne 19 Ne(n,p) 19 F 

with the result that about 70% of the 19 F produced via spallation of neutrons is 
then immediate destroyed, primarily by the (n, a) reaction above. In the next 
10~ 6 sec the coproduced protons are also processed 

15 N(p,a) 12 C 19 F(p,«) 16 23 Na(p,a) 20 Ne 

with the latter two reactions competing as the primary proton poisons. This 
makes an important prediction: stars with high Na abundances should make 
more F, as the 23 Na acts as a proton poison to preserve the produced F. 

• Finally, there is one other destruction mechanism, the heating associated with 
the passage of the shock wave. It turns out the the F produced prior to shock 
wave passage can survive if it is in the outside half of the Ne shell. The reaction 

19 F( 7 ,a) 15 N 

destroys F for peak explosion temperatures exceeding 1.7 • 10 9 K. Such a tempera- 
ture is produced at the inner edge of the Ne shell by the shock wave heating, but 
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not at the outer edge. 

If all of this physics in handled is a careful network code that includes the 
shock wave heating and F production both before and after shock wave passage, 
the following are the results: 

[ i9 F/ 2o Ne]/[ i9 F/ 2o Ne]0 Thcavy ^(MeV) 
0.14 4 
0.6 6 
1.2 8 
1.1 10 
1.1 12 

One sees that the attribution of F to the neutrino process argues that the heavy 
flavor v temperature must be greater than 6 MeV, a result theory favors. One 
also sees that F cannot be overproduced by this mechanism: although the instan- 
taneous production of F continues to grow rapidly with the neutrino temperature, 
too much F results in its destruction through the (p, a) reaction, given a solar 
abundance of the competing proton poison 23 Na. Indeed, this illustrates an odd 
quirk: although in most cases the neutrino process is a primary mechanism, one 
needs 23 Na present to produce significant F. Thus in this case the neutrino process 
is a secondary mechanism. 

While there are other significant neutrino process products ( 7 Li, 138 La, 
180 Ta, 15 N ...), the most important product is n B, produced by spallation off 
carbon. A calculation by Timmes et al. [18] found that the combination of the 
neutrino process, cosmic ray spallation and big-bang nucleosythesis together can 
explain the evolution of the light elements. The neutrino process, which produces 
a great deal of n B but relatively little 10 B, combines with the cosmic ray spalla- 
tion mechanism to yield the observed isotope ratio. Again, one prediction of this 
picture is that early stars should be n B rich, as the neutrino process is primary 
and operates early in our galaxy's history; the cosmic ray production of 10 B is 
more recent. There is hope that HST studies will soon be able to descriminate 
between 10 B and n B: as yet this has not been done. 

6. The r-process 

Beyond the iron peak nuclear Coulomb barriers become so high that charged 
particle reactions become ineffective, leaving neutron capture as the mechanism 
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responsible for producing the heaviest nuclei. If the neutron abundance is modest, 
this capture occurs in such a way that each newly synthesized nucleus has the 
opportunity to f3 decay, if it is energetically favorable to do so. Thus weak equi- 
librium is maintained within the nucleus, so that synthesis is along the path of 
stable nuclei. This is called the s- or slow-process. However a plot of the s-process 
in the (N,Z) plane reveals that this path misses many stable, neutron-rich nuclei 
that are known to exist in nature. This suggests that another mechanism is at 
work, too. Furthermore, the abundance peaks found in nature near masses A ~ 
130 and A ~ 190, which mark the closed neutron shells where neutron capture 
rates and j3 decay rates are slower, each split into two subpeaks. One set of sub- 
peaks corresponds to the closed-neutron-shell numbers N ~ 82 and N ~ 126, and 
is clearly associated with the s-process. The other set is shifted to smaller N, ~ 
76 and ~ 116, respectively, and is suggestive of a much more explosive neutron 
capture environment where neutron capture can be rapid. 

This second process is the r- or rapid-process, characterized by: 

• The neutron capture is fast compared to (5 decay rates. 

• The equilibrium maintained within a nucleus is established by (n, 7) <-> (7,71): 
neutron capture fills up the available bound levels in the nucleus until this equi- 
librium sets in. The new Fermi level depends on the temperature and the relative 
71/7 abundance. 

• The nucleosynthesis rate is thus controlled by the (3 decay rate: each f3~ cap- 
ture coverting n — > p opens up a hole in the neutron Fermi sea, allowing another 
neutron to be captured. 

• The nucleosynthesis path is along exotic, neutron-rich nuclei that would be 
highly unstable under normal laboratory conditions. 

• As the nucleosynthesis rate is controlled by the (5 decay, mass will build up at 
nuclei where the (5 decay rates are slow. It follows, if the neutron flux is reason- 
able steady over time so that equilibrated mass flow is reached, that the resulting 
abundances should be inversely proportional to these f3 decay rates. 

Let's first explore the (n, 7) <-* (7, n) equilibrium condition, which requires 
that the rate for (n, 7) balances that for (7, n) for an average nucleus. So consider 
the formation cross section 

A + rw(A+l)+7 



This is an exothermic reaction, as the neutron drops into the nuclear well. Our 
averaged cross section, assuming a resonant reaction (the level density is high in 
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heavy nuclei) is (see any standard nuclear astrophysics text, such as Clayton [6]) 



2tt \ 3/2 r n r 



=li^J '^ E ' KT (12) 

where E ~ is the resonance energy, and the Ts are the indicated partial and 
total widths. Thus the rate per unit volume is 

r <->~^(^f^ (13) 

where N n and Na are the neutron and nuclear number densities and \x the reduced 
mass. This has to be compared to the (7, n) rate. 

The (7, n) reaction requires the photon number density in the gas. This is 
given by the Bose-Einstein distribution 

nr/ x 87T e 2 de ,„ A . 

N(e) = -tt^^u^ — - 14 
v ' c 3 h 3 e e / kT - 1 

The high-energy tail of the normalized distribution can thus be written 

1 

where in the last expression we have set K — c — 1. 



Now we need the resonant cross section in the (7, n) direction. For photons 
the wave number is proportional to the energy, so 



ii r r 



As the velocity is c =1, 



N 1 vr 2 A^ 7 io e 2 (e - E r ) 2 + (r/2) 



2 -' /kT de- ^ ^f" 2 (16) 



We evaluate this in the usual way for a sharp resonance, remember that the energy 
integral over just the denominator above (the sharply varying part) is 2n/T: 



^l^n Er/kT ^ 

~ iv 7 r 



So that the rate becomes 



r M ~2N A J-^e~ E ^T (17) 
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Equating the (n, 7) and (7,71) rates and taking Na ~ ^Va_i then yields 



2 //xc 2 fcT 



) 



3/2 



n 



(he) 3 [ 2tt 



e -E r /kT 



(18) 



where the /is and cs have been properly inserted to give the right dimensions. 
Now E r is esssentially the binding energy. So plugging in the conditions N n ~ 
3 x 10 23 /cm 3 and T 9 ~ 1, we find that the binding energy is ~ 2.4 MeV. Thus 
neutrons are bound by about 30 times kT, a value that is still small compared to 
a typical binding of 8 MeV for a normal nucleus. (In this calculation I calculated 
the neutron reduced mass assuming a nuclear target with A=150.) 

The above calculation fails to count spin states for the photons and nuclei 
and is thus not quite correct. But it makes the essential point: the r- process 
involves very exotic species largely unstudied in any terrestrial laboratory. It is 
good to bear this in mind, as in the following section we will discuss the responses 
of such nuclei to neutrinos. Such responses thus depend on the ability of theory 
to extrapolate responses from known nuclei to those quite unfamiliar. 

The path of the r-process is along neutron-rich nuclei, where the neutron 
Fermi sea is just ~ (2-3) MeV away from the neutron drip line (where no more 
bound neutron levels exist). After the r-process finishes (the neutron exposure 
ends) the nuclei decay back to the valley of stability by (3 decay. This can involve 
some neutron spallation (/3-delayed neutrons) that shift the mass number A to 
a lower value. But it certainly involves conversion of neutrons into protons, and 
that shifts the r-process peaks at N ~ 82 and 126 to a lower N, off course. This 
effect is clearly seen in the abundance distribution: the r-process peaks are shifted 
to lower N relative to the s-process peaks. This is the origin of the second set of 
"subpeaks" mentioned at the start of the section. 

It is believed that the r-process can proceed to very heavy nuclei (A ~ 270) 
where it is finally ended by /^-delayed and n-induced fission, which feeds matter 
back into the process at an A ~ A max /2. Thus there may be important cycling 
effects in the upper half of the r-process distribution. 

What is the site(s) of the r-process? This has been debated many years 
and still remains a controversial subject. 
• The r-process requires exceptionally explosive conditions 



p(n) ~ 10 20 



cm 



3 T ~ 10 9 K t ~ Is. 
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• Both primary and secondary sites proposed. Primary sites are those not re- 
quiring preexisting metals. Secondary sites are those where the neutron capture 
occurs on preexisting s-process seeds. 

• Suggested primary sites include the the neutronized atmosphere above the proto- 
neutron star in a Type II supernova, neutron-rich jets produced in supernova ex- 
plosions or in neutron star mergers, inhomogeneous big bangs, etc. 

• Secondary sites, where p(n) can be lower for successful synthesis, include the 
He and C zones in Type II supernovae, the red giant He flash, etc. 

The balance of evidence favors a primary site, so one requiring no preen- 
richment of heavy s-process metals. Among the evidence: 

1) HST studies of very-metal-poor halo stars: The most important evidence are 
the recent HST measurements of Sneden et al. [15] of very metal-poor stars 
([Fe/H] ~ -1.7 to -3.12) where an r-process distribution very much like that of 
our sun has been seen for Z > 56. Furthermore, in these stars the iron content 
is variable. This suggests that the "time resolution" inherent in these old stars is 
short compared to galactic mixing times (otherwise Fe would be more constant). 
The conclusion is that the r-process material in these stars is most likely from 
one or a few local supernovae. The fact that the distributions match the solar 
r-process (at least above charge 56) strongly suggests that there is some kind of 
unique site for the r-process: the solar r-process distribution did not come from 
averaging over many different kinds of r-process events. Clearly the fact that 
these old stars are enriched in r-process metals also strongly argues for a primary 
process: the r-process works quite well in an environment where there are few 
initial s-process metals. 

2) There are also fairly good theoretical arguments that a primary r-process oc- 
curring in a core-collapse supernova might be viable [20]. First, galactic chemical 
evolution studies indicate that the growth of r-process elements in the galaxy is 
consistent with low-mass Type II supernovae in rate and distribution. More con- 
vincing is the fact that modelers have shown that the conditions needed for an 
r-process (very high neutron densities, temperatures of 1-3 billion degrees) might 
be realized in a supernova. The site is the last material blown off the supernova, 
the material just above the mass cut. When this material is blown off the star ini- 
tially, it is a very hot neutron-rich, radiation-dominated gas containing neutrons 
and protons, but an excess of the neutrons. As it expands off the star and cools, 
the material first goes through a freezeout to a particles, a step that essentially 
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locks up all the protons in this way. Then the as interact through reactions like 

a + a + a ^ u C 
a + a + n Be 

to start forming heavier nuclei. Note, unlike the big bang, that the density is 
high enough to allow such three-body interactions to bridge the mass gaps at A 
= 5,8. The a capture continues up to heavy nuclei, to A ~ 80, in the network 
calculations. The result is a small number of "seed" nuclei, a large number of as, 
and excess neutrons. These neutrons preferentially capture on the heavy seeds 
to produce an r-process. Of course, what is necessary is to have ~ 100 excess 
neutrons per seed in order to successfully synthesize heavy mass nuclei. Some of 
the modelers find conditions where this almost happens. 

There are some very nice aspects of this site: the amount of matter ejected 
is about 10~ 5 — 10~ 6 solar masses, which is just about what is needed over the 
lifetime of the galaxy to give the integrated r-process metals we see, taking a 
reasonable supernova rate. But there are also a few problems, especially the fact 
that with calculated entropies in the nucleon soup above the proto-neutron star, 
neutron fractions appear to be too low to produce a successful A ~ 190 peak. 
There is some interesting recent work invoking neutrino oscillations to cure this 
problem: charge current reactions on free protons and neutrons determine the 
n/p ratio in the gas. Then, for example, an oscillation of the type v e — > f s teriie 
can alter this ratio, as it would turn off the z/ e s that destroy neutrons by charged- 
current reactions. Unfortunately, a full discussion of such possibilities would take 
us too far afield today. 

The nuclear physics of the r-process tells us that the synthesis occurs 
when the nucleon soup is in the temperature range of (3-1) -10 9 K, which, in 
the hot bubble r-process described above, corresponds to a freezeout radius of 
(600-100) km and a time ~ 10 seconds after core collapse. The neutrino fluence 
after freezeout (when the temperature has dropped below 10 9 K and the r-process 
stops) is then ~ (0.045-0.015) -10 51 ergs/(100km). Thus, after completion of the 
r-process, the newly synthesized material experiences an intense flux of neutrinos. 
This brings up the question of whether the neutrino flux could have any effect on 
the r-process. 
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7. Neutrinos and the r-process 

Rather than describe the exotic effects of neutrino oscillations on the r- 
process, mentioned briefly above, we will examine standard-model effects that 
are nevertheless quite interesting. The nuclear physics of this section - neutrino- 
induced neutron spallation reactions - is also relevant to recently proposed su- 
pernova neutrino observatories such as OMNIS and LAND. In contrast to our 
first discussion of the //-process in Section 5, it is apparent that neutrino effects 
could be much larger in the hot bubble r-process: the synthesis occurs much 
closer to the star than our Ne radius of 20,000 km: estimates are 600-1000 km. 
The r-process is completed in about 10 seconds (when the temperature drops to 
about one billion degrees), but the neutrino flux is still significant as the r-process 
freezes out. The net result is that the "post-processing" neutrino fluence - the 
fluence that can alter the nuclear distribution after the r-process is completed - 
is about 100 times larger than that responsible for fluorine production in the Ne 
zone. Recalling that 1/300 of the nuclei in the Ne zone interacted with neutrinos, 
and remembering that the relevant neutrino-nucleus cross sections scale as A, 
one quickly sees that the probability of a r-process nucleus interacting with the 
neutrino flux is approximately unity. 

Because the hydrodynamic conditions of the r-process are highly uncer- 
tain, one way to attack this problem is to work backward in time. We know the 
final r-process distribution (what nature gives us) and we can calculate neutrino- 
nucleus interactions relatively well. Thus from the observed r-process distribution 
(including neutrino postprocessing) we can work backward to find out what the 
r-process distribution looked like at the point of freezeout. In Figs. 5 and 6, the 
"real" r-process distribution - that produced at freezeout - is given by the dashed 
lines, while the solid lines show the effects of the neutrino postprocessing for a 
particular choice of fluence [12]. The nuclear physics input into these calculations 
is precisely that previously described: GT and first-forbidden cross sections, with 
the responses centered at excitation energies consistent with those found in ordi- 
nary, stable nuclei, taking into account the observed dependence on \N — Z\. 

One important aspect of the figures is that the mass shift is significant. 
This has to do with the fact that a 20 MeV excitation of a neutron-rich nucleus 
allows multiple neutrons ( ~ 5) to be emitted. (Remember we found that the 
binding energy of the last neutron in an r-process neutron-rich nuclei was about 
2-3 MeV under typical r-process conditions.) The second thing to notice is that 
the relative contribution of the neutrino process is particularly important in the 
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Fig. 5. Comparison of the r-process distribution that would result from the 
freezeout abundances near the A ~ 130 mass peak (dashed line) to that where 
the effects of neutrino postprocessing have been include (solid line) . The fluence 
has been fixed by assuming that the A = 124-126 abundances are entirely due 
to the //-process. 

"valleys" beneath the mass peaks: the reason is that the parents on the mass 
peak are abundant, and the valley daughters rare. In fact, it follows from this 
that the neutrino process effects can be dominant for precisely seven isotopes (Te, 
Re, etc.) lying in these valleys. Furthermore if an appropriate neutrino fluence 
is picked, these isotope abundances are produced perfectly (given the abundance 
errors). The fluences are 

N = 82 peak 0.031 • 10 51 ergs/(100km) 2 /flavor 

N = 126 peak 0.015 • 10 51 ergs/(100km) 2 /flavor 

values in fine agreement with those that would be found in a hot bubble r-process. 
So this is circumstantial but significant evidence that the material near the mass 
cut of a Type II supernova is the site of the r-process: there is a neutrino finger- 
print. 

In conclusion, I hope this whirlwind tour through the nuclear aspects of 
neutrino interactions in detectors and in stars has illustrated how nuclear, neu- 
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Fig. 6. As in Fig. 5, but for the A ~ 195 mass peak. The A = 183-187 
abundances are entirely attributed to the //-process. 

trino, and stellar physics is interconnected. With the HST and other great obser- 
vatories providing so much new information on the nuclear microphysics governing 
the universe, it is likely that astrophysics will continue to provide important chal- 
lenges to nuclear physicists. Likewise, examples like the "neutrino fingerprint" on 
the r-process illustrate how an understanding of nuclear physics can help astro- 
physicists settle issues like the site of the r-process. 

The work was supported in part by the US Department of Energy. 
8. References 

1. E. G. Adelberger and W. C. Haxton, 1987, Phys. Rev. C 36, 879. 

2. Austin, S. M., Anantaraman, N., and Love, W. G., 1994, Phys. Rev. Lett. 73, 
30; Watson, J. W. et al., 1985, Phys. Rev. Lett. 55, 1369. 

3. Bahcall, J. N., 1989, Neutrino Astrophysics (Cambridge University Press, 
Cambridge). 

4. Bahcall, J. N. and Barnes, C. A., 1964, Phys. Lett 12, 48. 

5. Bethe, H. and Wilson, J. R., 1985, Ap. J. 295, 14. 

6. Clayton, D. D., 1968, Principles of Stellar Evolution and Nucleosynthesis 



30 



(McGraw-Hill, New York). 

7. Cleveland, B. T. et al, 1998, Ap. J. 496, 505. 

8. Cooperstein, J., Bethe, H. A., and Brown, G. E., 1984, Nucl. Phys. 429, 527. 

9. Fuller, G. M. and Meyer, B. S., 1995, Ap. J. 453, 792. 

10. Garcia, A. et al, 1991, Phys. Rev. Lett. 67, 3658 and 1990, Phys. Rev. C 
42, 765; Trindler, W. et al., 1995, Phys. Lett. B 349, 267. 

11. Goldhaber, M. and Teller, E., 1948, Phys. Rev. 74, 1046; Donnelly, T. W., 
Dubach, J., and Haxton, W. C, 1975, Nucl. Phys. A 251, 353. 

12. Haxton, W. C, Langanke, K., Qian, Y.-Z., and Vogel, P., 1997, Phys. Rev. 
Lett. 78, 2694 and Phys. Rev. C 55, 1532. 

13. Mezzacappa, A. et al., 1998, Ap. J. 495, 911; Janka, H.-Th. and Muller, E., 
1996, Astron. Astrophys. 306, 167; Burrows, A., Hayes, S., and Fryxell, B. 
A., 1995, Ap. J. 450, 830. 

14. Nakayama, K., Galeao, A. P., and Krmpotic, K., 1982, Phys. Lett. B 114, 
217; D. J. Horen et al., 1980, Phys. Lett. B 95, 27. 

15. Pfeiffer, B. et al., 1998, |astro-ph/ 98 12414 . 



16. Poskanzer, A. M., McPherson, R., Esterlund, R. A., and Reeder, P. L., 1966, 
Phys. Rev 152, 995; Sextro, R. G., Gough, R. A., and Cerny, J, 1974, Nucl. 
Phys. A 234, 130. 

17. Rapaport, J. et al., 1981, Phys. Rev. Lett. 47, 1518. 

18. Timmes, F. X., Woosley, S. E., and Weaver, T. A., 1995, Ap. J. Suppl. 98, 
617. 

19. Walecka, J. D., 1975, in Muon Physics, ed. V. W. Hughes and C. S. Wu 
(Academic Press, New York), vol. 2, p. 113. 

20. Woosley, S. E. et al., 1994, Ap. J. 433, 229; Woosley, S. E. and Hoffman, R. 
D., 1992, Ap. J. 395, 202. 

21. Woosley, S. E. and Haxton, W. C, 1988, Nature 334, 45; Woosley, S. E., 
Hartmann, D. H., Hoffman, R. D., and Haxton, W. C, 1990, Ap. J. 356, 272. 



